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Abstract 

There can exist a parallel 'mirror' world which has the same particle physics as the 
observable world and couples the latter only gravitationally. The nucleosynthesis bounds 
demand that the mirror sector should have a smaller temperature than the ordinary one. 
By this reason its evolution should be substantially deviated from the standard cosmology 
as far as the crucial epochs like baryogenesis, nucleosynthesis etc. are concerned. Starting 
from an inflationary scenario which could explain the different initial temperatures of the two 
sectors, we study the time history of the early mirror universe. In particular, we show that 
in the context of the GUT or electroweak baryogenesis scenarios, the baryon asymmetry in 
the mirror world should be larger than in the observable one and in fact the mirror baryons 
could provide the dominant dark matter component of the universe. In addition, analyzing 
the nucleosynthesis epoch, we show that the mirror helium abundance should be much larger 
than that of ordinary helium. The implications of the mirror baryons representing a kind 
of self-interacting dark matter for the large scale structure formation, the CMB anysotropy, 
the galactic halo structures, microlensing, etc. are briefly discussed. 
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1 Introduction 



The old idea that there can exist a hidden mirror sector of particles and interactions which is 
the exact duplicate of our visible world has attracted a significant interest over last years, in 
particular being motivated by the problems of neutrino physics [||, |j| , gravitational microlensing 
j|, ||], gamma ray bursts H, ultra- high energy cosmic rays (7j, flavour and CP violation M, 
etc. The basic concept is to have a theory given by the product G x G' of two identical gauge 
factors with the identical particle contents, which could naturally emerge e.g. in the context of 
E% x E' 8 superstring. Two sectors communicate through gravity and perhaps also via some other 
messengers. A discrete symmetry P(G <-> G') interchanging corresponding fields of G and G' , so 
called mirror parity, implies that both particle sectors are described by the same LagrangiansJJ 

In particular, one can consider a minimal symmetry Gsm x G' sm where Gsm = SU{3) x 
SU(2) x U(l) stands for the standard model of observable particles: three families of quarks and 
leptons qi, u\, d\\ k, (i = 1,2,3) and the Higgs doublet (ft, while G' SM = [SU(3) x SU(2) x 
U(l)]' is its mirror gauge counterpart with analogous particle content: fermions q' iri u'f, df; l[, e'f 
and the Higgs (/>'. (From now on all fields and quantities of the mirror (M) sector will have an 
apex to distinguish from the ones belonging to the observable or ordinary (O) world.) The 
mirror parity implies that all coupling constants (gauge, Yukawa, Higgs) have the same pattern 
in both sectors and thus their microphysics is the same.0 

One could naively think that due to mirror parity the O- and M- particles should have 
the same cosmological densities, which would be in the immediate conflict with the Big Bang 
nucleosynthesis (BBN) bounds on the effective number of extra light neutrinos, AN U < 1 [11]: 
the mirror photons, electrons and neutrinos would give a contribution to the Hubble expansion 
rate equivalent to AN U ~ 6.14. Therefore, the M-particle density in the early universe should 
be appropriately reduced. This situation is plausible if two following conditions are satisfied: 

A. At the initial moment two systems are born with different densities. In particular, the 
inflationary reheating temperature in the M-sector should be lower than in the visible one, 
T' R < Tr, which can be achieved in certain models [7|, [Tcfl . 

B. The M- and O-particles interact very weakly, so that two systems do not come into 
the thermal equilibrium with each other in the early universe. This condition is automatically 
fulfilled if two worlds communicate only via the gravity. More generally, there could be other 
messengers like superheavy gauge singlet fields or light singlets of the moduli type. In either 
case, they should mediate the effective couplings between the O- and M- particles suppressed 
by a large mass factor M ~ Mp or so. 

If two sectors have different reheating temperatures, and if they do not come into the thermal 
contact at later stages, then during the universe expansion they evolve independently and ap- 
proach the BBN epoch with different temperatures. Namely, the BBN bound AN U < 1 implies 
that T'/T < 0.64. 

In this paper we study the comparative time history of two sectors in the early universe. 
We show that due to the temperature difference, in the mirror sector all key epochs as are 
the baryogenesis, nucleosynthesis, etc. proceed at somewhat different conditions than in the 
observable universe. The paper is organized as follows. In sect. 2 we describe an inflationary 
model which could provide different reheating temperatures between two sectors. The sect. 3 we 
show that the baryon asymmetry in the M-world generically should be larger than that in the O- 
sector either in the GUT or electroweak baryogenesis scenarios. Moreover, it is pretty plausible 
that M-baryons provide a significant fraction of the dark matter of the universe, i.e. £l' B > £Ib- In 

x In the brane world picture, the M-sector can be the same O-world realized on a parallel brane, G' = G ||. 
2 The mirror parity could be spontaneously broken and the weak interaction scales '•((/)) = V and {</>') = v' could 
be different, which leads to somewhat different particle physics in the mirror sector Q 
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sect. 4 we study the chemical composition of the M-sector and find that it should be dominantly 
a helium world - the primordial abundance of the mirror 4 He should be 2 — 3 times larger than 
of the ordinary one. In sect. 5 we briefly address the problem of the cosmological structure 
formation in the presence of M-baryons as a dominant dark matter component. Namely, M- 
baryons being a sort of self-interacting dark matter could provide interesting signatures on 
the CMB anisotropy, the large scale structure of the universe, the form of the galactic halos, 
microlensing, etc. Finally, in sect. 6 we summarize our findings. 



2 Inflation and post-inflation 

An attractive realization of the inflationary paradigm is provided by supersymmetric models of 
hybrid inflation. The symplest model is based on the superpotential W = \S(& 2 — fi 2 ) containing 
the inflaton field S and the additional 'orthogonal' field <3?, where A is order 1 coupling constant 



and /i is a dimensional parameter of the order of the GUT scale 13]. The supersymmetric 
vacuum is located at S = 0, $ = fj,, while for the field values <J> = 0, S > /i the tree level 
potential has a flat valley with an energy density V = A 2 /i 4 . Since the supersymmetry is broken 
by the non- vanishing F-term, F$ = A/i 2 , the flat direction is lifted by radiative corrections and 
the potential of S gets a slope which is appropriate for the slow roll conditions. The COBE 
results on the CMB anisotropy imply that y 1 / 4 ~ e 1 / 4 x 7 • 10 16 GeV, where e <C 1 is a slow-roll 
parameter. 

The genetic problem of the hybrid inflation models concerns the choice of the initial condi- 
tions [14, |i~q]. Namely, at the end of the Planck era the scalar S should have an initial value 



~ Mp while $ must be zero with the high accuracy over a region much larger than the initial 
horizon size ~ Mp 1 . In other words, the initial field configuration should be located just on the 
bottom of the S valley, and so the energy density to start with should be ~ /i 4 , at least eight 
orders of magnitude smaller than the natural energy density Mp at the Planck era. If the <E> 
field would have initial amplitude ~ Mp, then its oscillation should be damped within typical 
time ~ Mp , otherwise it induces too a big curvature for the inflaton S and violate the slow-roll 
conditions, thus preventing the onset of the inflation. 

Such a Fine Tuning can be avoided by a possible preheating before the onset of the inflation: 
the oscillating field could decay into some light particles |]l5[ . In this view, one can consider 
a model based on superpotential (ll| : 

^preheat = \S(<S> 2 - /J 2 ) + 5 $(tf 2 + *' 2 ) (1) 

Now all fields can have initial values order Mp, and thus the initial energy density can be ~ Mp. 
At first instants, due to oscillation of <£, the system behaves as matter dominated universe. Fast 
damping of $ is a pretext of inflationary stage which allows the inflaton energy density ~ fi 4 
to dominate and then S can slowly roll to the origin. This function is carried by the second 
term in (|l]) — the oscillating orthogonal field $ fastly decays into ^ and *S>' particles which have 
practically no contact to inflaton S. In addition, with vanishing also effective mass terms of 
^'s disappear and the latter fields start to behave as massless - they stop oscillating and freeze. 
In general, oscillations and ^' freeze at different amplitudes, typically ~ /i, at which they are 
catched by the moment when their mass g& drops below the Hubble parameter (see Fig. 3 of 



ref. jl5|]). When slow roll ends up, all fields start to oscillate around their vacuum values, 5 = 0, 
<I> = ji, ^' = 0, and reheat the universe. 

Let us assume now that reheating occurs due to superpotential terms: 

Wreheat = hVfafc + KV ' M'2 ( 2 ) 
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where (p\ t 2 and <j)'± 2 are the Higgs doublet superfields respectively for the O- and M-sectors. 
Then different magnitudes of \P and W at the end of slow-roll phase should reflect into difference 
of reheating temperatures Tr and T' R in two systems, simply because of the energy difference 
stored into and oscillations which decay respectively into O- and M-Higgses. More detailed 
analysis of this mechanism will be presented elsewhere.^ 

We have also to make sure that after reheating two sectors do not come into the thermal 
equilibrium to each other, or in other words, that interactions between O- and M-particles are 
properly suppressed. In a supersymmetric theory this condition is fulfilled in a rather natural 
manner. In particular, the mixed terms of the lowest dimension in the superpotential are: 

+ (3) 

which are suppressed by a large mass factor M ~ Mp or so, and thus are safe. The same 
holds true for soft supersymmetry breaking F- and D-terms like [jj((/>i^2)( < P'i ( P2)]f, etc., where 
z = msO 2 is the supersymmetry breaking spurion with ms ~ 1 TeV. As for the the kynetic 
mixing term F^ v F' between the field-strength tensors of the gauge factors U(l) and U(l)', it 
can be forbidden by embedding Gsm x G' sm in the grand unified group like SU(5) x SU(5)'. 

Once the O- and M-systems are decoupled already after reheating, at later times t they will 
have different temperatures T(t) and T'(t), and so different energy and entropy densities: 

Pit) = ^9*(T)T\ p'(t) = ^(T')T' 4 , (4) 

s(t) = 2 -^9s{T)T\ s'(t) = ^(T')T' 3 . (5) 

The factors g s and g[, g' s accounting for the effective number of the degrees of freedom in 
two systems can in general be different from each other. During the universe expansion, the two 
sectors evolve with separately conserved entropies. Therefore, the ratio x = (s'/s) 1 / 3 is time 
independent,^] while the ratio of the temperatures in two sectors is simply given by: 



T(t) ~ 



9s{T)f'\ (6) 



SfsF') 



The Hubble expansion rate is determined by the total energy density p = p + p', H 
\J (8tt/3)Gnp. Therefore, at a given time t in a radiation dominated epoch we have 



T 
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H(t) = - = 1.66V<7*(T)— = 1.66^')^ (7) 
in terms of O- and M-temperatures T(t) and T'{t), where 

UT) = 5 ,(T)(1 + ax% gi(T') = g'^T') U + ■ (8) 

Here the factor a(T, T') = [g'*i T ') 1 9*i T )\ ' bs( T )/Ss( T ')] 4/3 takes into account that for T ^ T 
the relativistic particle contents of the two worlds can be different. However, except for very 
small values of x, we have a ~ I. So hereafter we always take g*(T) = 5*(T)(I + x 4 ) and 
g'*iT') = g£(T')(l + x -4 ). In particular, in the modern universe we have a(To,Tg) = 1, <7 s (Tb) = 

3 For other scenarios of the asymmetric reheating, see ^,0, pTj| . 

4 We assume that expansion goes adiabatically in both sectors and neglect the additional entropy production 
due to the possible weakly first order electroweak or QCD phase transitions. 
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g' s (To) = 3.91, and x = Tq/Tq, where Tq,Tq are the present temperatures of the 0- and M- relic 
photons .0 In fact, x is the only free parameter in our model and it is constrained by the BBN 
bounds. 

The observed abundances of light elements are in good agreement with the standard nucle- 
osynthesis predictions, when at T ~ 1 MeV we have g* = 10.75 as it is saturated by photons 
7, electrons e and three neutrino species v e ,u,T- The contribution of mirror particles (7', e' and 
^e,ur) would change it to g* = g*(l + x 4 ). Deviations from g* = 10.75 are usually parametrized 
in terms of the effective number of extra neutrino species, Ag = g* — 10.75 = 1.75AJVj,. Thus 
we have: 

AiV„ = 6.14 • x 4 . (9) 



In view of the present observational situation, a reliable bound is AN U < 1 [11], which translates 
as x < 0.64. This limit very weakly depends on AN U . E.g. AN V < 1.5 implies x < 0.70. 

As far as x 4 <C 1, in a relativistic epoch the Hubble expansion rate (|7|) is dominated by the O- 
matter density and the presence of M-sector practically does not affect the standard cosmology 
of the early ordinary universe. However, even if the two sectors have the same microphysics, 
the cosmology of the early mirror world can be very different from the standard one as far as 
the crucial epochs like baryogenesis, nuclesosynthesis, etc. are concerned. Any of these epochs 
is related to an instant when the rate of the relevant particle process T(T), which is generically 
a function of the temperature, becomes equal to the Hubble expansion rate H(T). Obviously, 
in the M-sector these events take place earlier than in the O-sector, and as a rule, the relevant 
processes in the former freeze out at larger temperatures than in the latter. 

In the matter domination epoch the situation becomes different. In particular, we know 
that ordinary baryons can provide only a minor fraction of the present cosmological density, 
0,B = 0.01 — 0.06, whereas the observational data indicate the presence of dark matter which 
can amount for Q, m ~ 0.2 — 1. It is interesting to question whether the missing matter density 
of the universe could be due to mirror baryons? In the next section we show that this situation 
can emerge in a pretty natural manner. 



3 Baryogenesis and mirror baryon density 

It is well known that a non-zero baryon asymmetry (BA) can be produced in the initially 
baryon symmetric universe if three following conditions are fulfilled: B-violation, CP violation 
and departure from the thermal equilibrium. Generally speaking, the baryogenesis scenarios can 
be divided in two categories in which the out of equilibrium conditions are provided (a) by the 
universe expansion itself, or (b) by fast phase transition and bubble nucleation. In particular, 
the latter concerns the electroweak baryogenesis schemes, while the former is typical for a GUT 
type baryogenesis or leptogenesis. 

At present it is not possible to say which of the known mechanisms is responsible for the 
observed BA. We only know that the baryon to photon number density ratio r] = ns/n^ is 
restricted by the BBN constraints to the range 77 = (2 - 6) x 10~ 10 . It is most likely that the 
BA in the M- world rj' = n' B /nL is produced by same mechanism and moreover, the rates of the 
B and CP violation processes are parametrically the same in both cases. However, the out of 
equilibrium conditions should be different since at relevant temperatures the universe expansion 

5 The frozen ratio of the neutrino and photon temperatures in the M-sector r' = T^ /Tq depends on the 1/ 
decoupling temperature from the mirror plasma which scales approximatively as T' D ~ x~ 2 ^ 3 Td, where To = 2 — 3 
MeV is the decoupling temperature of the usual neutrinos. Therefore, unless x < 10 -3 , r' has a standard value 
ro = Tuo/To = (4/ll) 1/3 . For x < 1(T 3 , T' D becomes larger than the QCD scale A ~ 200 MeV, so that due to 
the mirror gluons and light quarks u',d',s' contribution we would obtain r' = (4/53) 1 / 3 , g' s (To) = 2.39 and so 
Tq/To — 1.2x. However, in the following such small values of x are not of our interest. 
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is faster for the M-sector. Below we show that by this reason rf typically emerges larger than 77 
for either type (a) or (b) scenarios. 

The M-baryons can be of the cosmological relevance if £l' B exceeds Qb = 3.67 x 10 7 rjh~ 2 = 
0.01 — 0.06, whereas Q' B > 1 would overdose the universe. So we are interested in a situation 
when the ratio f3 = Q' b /Qb falls in the range from 1 to about 100. Since = x 3 n 7 , we obtain 
(3 = x 3 r]' /rj. Therefore, rf > rj does not apriori mean that (5 > 1, and in fact there is a lower 
limit x > 10~ 2 or so for the relevant parameter space. Indeed, it arises from x 3 = f3r]/rj' by 
recalling that r) ~ 10~ 9 , while rf can be taken at most ~ 10 -3 , the biggest value which can be 
principally realized in any baryogenesis scheme under the realistic assumptions. 



3.1 GUT Baryogenesis 

The GUT baryogenesis mechanism typically based on a superheavy boson X undergoing the B- 
and CP-violating decays into quarks and leptons. The following reaction rates are of relevance: 
Decay: To ~ axMx for T J$ Mx or ~ axM x /T for T ^ Mx, where ax is the coupling 
strength of X to fermions and Mx is its mass; 

Inverse decay: T ID ~ T D (M X /T) 3 / 2 exp(-M x /T) for T < M x or T ID ~ T D for T > M x ; 
The X boson mediated 2^2 processes: T$ ~ rixcr ~ Aa 2 x T 5 /(M x + T 2 ) 2 , where the factor A 
amounts for the possible reaction channels. 

The final BA depends on a temperature at which X bosons go out from equilibrium. One can 



introduce a parameter which measures the effectiveness of the decay at the epoch T ~ Mx UW - 

— 1/2 

k = (Td/2H)t=m x = 0.3g* (axMpijMx)- The larger is k the longer equilibrium is main- 
tained and the freeze-out abundance of X boson becomes smaller. Hence, the resulting baryon 
number to entropy ratio, B = tib/s — 0.14r] is a decreasing function of k. It is approximately 
given as B ~ jj-F(k, k c ), where e is the CP violating factor and 

r 1 if k < 1 

F(k,k c ) = I 0.3k- 1 (log k)-° s ifl<k<k c (10) 

I VAa x ke-l( Aax ^ 1/4 if k > k c 

Here k c is a critical value defined by equation k c (logk c )~ 2A = 300/ (Aax)- It distinguishes 
between the regimes k < k c , in which inverse decay is relevant and k > k c , in which instead 
2 «-> 2 processes are the dominant reason for baryon damping. 

In a general context, without referring to a particular model, it is difficult to decide which 
range of parameters k and k c can be relevant for baryogenesis. One can impose only the most 
reasonable constraints g s (T = Mx) > 100 and e < 10 -2 , and thus e/g s < 10~ 4 or so. For a 
given mechanism responsible for the observed baryon asymmetry B ~ 10~ 10 , this translates into 
a lower bound F(k, k c ) > 10~ 6 . 

The presence of the mirror sector practically does not alter the ordinary baryogenesis. The 
effective particle number is g*(T) = g*(T)(l + x 4 ) and thus the contribution of M-particles to 
the Hubble constant at T ~ Mx is suppressed by a factor x 4 . 

In the mirror sector everything should occur in a similar way, apart from the fact that now 
at T' ~ Mx the Hubble constant is not dominated by the mirror species but by ordinary ones: 
g'*(T') ~ <^(T")(1 + x~ 4 ). As a consequence, we have k! = (T' D /2H)\ T / =Mx = kx 2 . Since the 
value of k c is the same in the two sectors, the mirror baryon asymmetry can be simply obtained 



by replacing k — > k' = kx in eq.10, i.e. B' = n' B /s' ~ (e/g' s )F(k' = kx ,k c ). Since F is a 



decreasing function of k, then for x < 1 we have F(kx 2 ,k c ) > F(k,k c ) and thus we conclude 
that the mirror world always gets a larger BA than the visible one, B' > B. 

However, this does not apriori mean that Vt' B is always larger than Qb- Since the entropy 
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Figure 1: Panel A. The combination x 3 F(kx 2 ,k c ) as a function of k for k c = 10 4 and x = 0.6, 0.1, 0.01 
(dash). The solid curve corresponding to x = 1 in fact measures the possible BA in the ordinary world, 
F(k, k c ) = (g s /e)B(k). Panel B. The curves confining the parameter region in which f3 = £1' b /Qb varies 
from 1 to 100, for x — 0.6 (solid) and for x = 0.01 (dash). The parameter area above thick solid curve 
corresponds to F(k, k c ) < 10~ 6 and it is excluded by the observable value of r\. 



densities are related as s'/s = x , for the ratio f3 = VL' b /£Ib we have: 

The behaviour of the factor x^F(kx 2 , k c ) as a function of k for different values of the parameter 
x is given in the Fig. 1A. Clearly, in order to have £l' B > Qb the function F(k,k c ) have to 
decrease faster than k~ 3 ^ 2 between k' = kx 2 and k. Closer inspection of the function [H] reveals 
that the M-baryons can be overproduced only if k is order k c or larger. In other words, the 
relevant interactions in the observable sector maintain equilibrium longer than in the mirror 
one, and thus ordinary BA can be suppressed by an exponential Boltzmann factor while the 
mirror BA could be produced still in non-exponential regime k' < k c . 

In Fig. IB we show the parameter region in which (3 = Q! b /Q,b falls in the range 1 — 100, 
in confront to the parameter area excluded by condition F(k,k c ) > 10~ 6 . We see that for 
x = 0.6 there is an allowed parameter space in which (3 can reach values up to 10, but (3 = 100 
is excluded. For a limiting case x = 10~ 2 , as it was expected, the parameter space for f3 > 1 
becomes incompatible with F(k,k c ) > 10~ 6 . For intermediate values of x, say x ~ 0.1 — 0.3, 
also the values (3 ~ 100 can be compatible. 

The above considerations can be applied also in the context of leptogenesis. One should 
remark, however, that potentially both the GUT baryogenesis or leptogenesis scenarios are in 
conflict with the supersymmetric inflation scenarios, because of the upper limit on the reheating 
temperatures about Tr < 10 9 GeV from the thermal production of gravitinos [17]. Moreover, it 



was shown recently that the non-thermal gravitino production can impose much stronger limits, 
Tr < 10 5 GeV or so p|. This probl em can be fully avoided m the electroweak baryogenesis 



scenario, which instead is definitely favoured by the supersymmetry. 
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3.2 Electroweak Baryogenesis 

The electroweak (EW) baryogenesis mechanism is based on the anomalous B-violating processes 
induced by the sphalerons which are quit rapid at high temperatures, but become much slower 
when temperature drops below 100 GeV. A succesfull scenario needs the first order EW phase 
transition and sufficient amount of CP violation, which conditions can be satisfied in the frames 
of the supersymmetric standard model, for certain parameter ranges [19]. 



The characteristic temperature scales of the electroweak phase transition are fixed entirely 
by the form of the finite temperature effective potential: 

V((j), T) = D(T 2 - T$)<j) 2 - ETcj) 3 + \ T (j) 4 (12) 

where all parameters can be expressed in terms of the fundamental couplings in the Lagrangian. 
For large temperatures, T S> 100 GeV, the electroweak symmetry is restored and V(4>, T) has 
a minimum at (j) = 0. With the universe expansion the temperature drops, approaching the 
specific values which define the sequence of the phase transition. These are all in the 100 GeV 
range and ordered as T\ > T c > T b > Tq. Namely, below T = T\ the potential gets a second 
local minimum 4> + (T). At the critical temperature T = T c the latter becomes degenerate with 
the false vacuum cj) = 0. At temperatures T < T c to the true vacuum state eft = 4>+(T) becomes 
energetically favoured, and transition to this state can occur via thermal quantum tunneling, 
through the nucleation of the bubbles which then expand fastly, percolate and finally fill the 
whole space within a horizon by the true vacuum. 

The bubble production starts when the free energy barrier separating the two minima be- 
comes small enough. The bubble nucleation temperature Tf, is defined as a temperature at which 
the probability for a single bubble to be nucleated within a horizon volume becomes order one: 

P(T b )=u(— ^- [i--±) e - n ~i (13) 



\H(T C )J V T t 

where F C (T) is the free energy and uj is an order one coefficient In particular, in the limit 
of thin wall approximation we have: 



F C (T) 64vr E (T c -T { 



81 (2A) 3 / 2 V T c 



(14) 



The condition (13) results into large values of F C (T{,)/T{„ typically order 10 2 . 

Once the bubble nucleation and expansion rate is larger than the Hubble parameter, the 
out of equilibrium condition for anomalous B-violating processes is provided by the fast phase 
transition itself. The BA can be produced inside the bubbles due to the CP violation since the 
quarks and antiquarks have different reflection coefficients on the bubble wall. The baryogenesis 
rate is completely independent from the universe expansion and it occurs practically in one 
instant as compared to the cosmological time scale of this epoch. 

As for the mirror sector, it is described by the same Lagrangian as the ordinary one and so 
the effective thermal potential of the mirror Higgs V(<p',T') has the same form as (|l2|). Then 
the temperature scales which are defined entirely by the form of the effective potential should 
be exactly the same for O- and M-sectors. Namely, T[ = T\ and T' c = T c . 

The equation ( |i~3| ) is the same for the M-sector apart of the fact that the corresponding 
Hubble constant is different: H(T' = T c ) = x~ 2 H(T = T c ). Therefore, we obtain: 

^7^ = ^^ + 81ogx (15) 
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Figure 2: The contours of /3 = Q' b /D,b m the plane of the parameters x and D, corresponding to 
(3 = 1, 10 and 100 from top to bottom. 



which in turn tells that the bubble nucleation temperatures in two sectors are practically equal: 
T' h = T b (l + 0.01 log x). (Clearly, the phase transition in M-sector occurs at earlier time than 
in O-sector: t' b ~ x~ 2 t b .) The reason is that between the temperature scales T c and To the free 
energy F C (T) is a rapidly changing function but the change in temperature itself is insignificant. 
Hence, we expect that the initial BA's produced right at the phase transition to be the same in 
O- and M-sectors: B(T = T b ) = B'{T' = T b ). 

However, the instantly produced baryon number can be still washed out by the sphaleron in- 
teractions. The anomalous B-violation rate T(T) ~ exp[— F(T)/T], where F(T) is the sphaleron 
free energy at finite T, may be large enough inside the bubble as far as the temperature is large. 
But it quickly falls as the temperature decreases, and baryon number freezes out as soon as T(t) 
drops below H(t) (^). The wash-out equation dB/dt = —T(t)B can be rewritten as 

dB F(T) , , 

~b = lrF dT < la > 

and integrated. Then the final BA in the O- and M-sectors can be expressed respectively as 

B = B(T b )D^ 1+x ^' 1/2 , B' = B(T b )D x ^ 1+x ^~ 1/2 , (17) 
where D < 1 is the baryon number depletion factor: 



D = exp 



-1/2 f b J 

-0.6^ 1 Mpi / dT- 



T 3 



(18) 



and g* ~ O(100) in the super symmetric standard model. Thus we always have B' > B, while 
the M-baryon mass density relative to O-baryons reads: 



q B 1 „, % 1 — x 2 



P(x) = ^- = x 3 ^- = x*D- K W, K{x) = ■ (19) 



n B b v ' y/fTx 4 

Lacking a precise theory for non-perturbative sphaleron transitions in the broken phase, the 
exact value of D one cannot be calculated even in the context of concrete models. If D ~ 1, 
the wash-out is ineffective and practically all BA produced right at the bubble nucleation is 
conserved. In this case £l' B should smaller than However, if D is enough small, one can 
achieve sufficiently large £l' B . The contourplot for the parameters x and D for which (3 falls 
in the range 1 — 100 is given in Fig. 2. For small x we have essentially j3 ~ x^D~ l and thus 
100 > /3 > 1 requires a depletion factor in the interval D = (10~ 2 — l)x 3 . Once again, for 
x ~ 10~ 2 one needs the marginal values D ~ 10 -8 — 10 -6 below which the observable BA 
B ~ 10 -10 cannot be produced at all. 
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4 Primordial Nucleosynthesis and mirror helium abundance 



The time scales relevant for standard BBN are denned by the "freeze-out" temperature of weak 
interactions Ty/ ~ 0.8 MeV (tw ~ 1 s) and by the "deuterium bottleneck" temperature Tjv — 
0.07 MeV (tjv ~ 200 s) fl6|| . When T > T\y, weak interactions transform neutrons into protons 
and viceversa and keep them in chemical equilibrium. The neutron abundance X n = n n /riB, 
defined as the ratio of neutron to baryon densities, is given by X n {T) = [1 + exp(Am/T)] _1 , 
where Am ~ 1.29 MeV is the neutron-proton mass difference. For T < Tyy the weak reaction 
rate T\y ~ G^T 5 drops below the Hubble expansion rate H(T) ~ 5.5T 2 /M Ph the neutron 
abundance freezes out at the equilibrium value X n (Tw) and it then evolves only due to the 
neutron decay: X n (t) = X n {Tw) exp(— i/r), where r = 886.7 s is the neutron lifetime. 

At temperatures T > T/v, the process p + n <-> d + 7 is faster than the universe expansion, 
and free nucleons and deuterium are in chemical equilibrium. The light element nucleosynthesis 
essentially begins when the system cools down to the temperature 

^ - h (,) + 1 I W T») - a0 ' M ' T ' < 2 °> 

where = 2.22 MeV is the deuterium binding energy, and tun is the nucleon mass. Below this 
temperature the deuterium abundance starts to grow which in turn allows to produce also the 
heavier nuclei. Nearly all neutrons present at this time are finally captured in He nuclei, due 
to the large binding energy of the latter. Thus, the primordial 4 He mass fraction is: 

n^2x„(w= ,f p( :,'i; ) ^ o^. (2D 

1 + exp(Am/Tyi/) 

As we have already discussed, the presence of the mirror sector with a temperature T' <C T 
has practically no impact the standard BBN. In fact, the limit x < 0.64 has been set by 
uncertainties of the present observational situation. In the mirror sector nucleosynthesis proceeds 
along the same lines. However, the impact of the O-world for the mirror BBN is dramatic. 

For any given temperature T", now we have H(T') ~ 5.5(1 + x~ 2 )T' 2 /Mpi for the Hubble 
expansion rate. Therefore, the freeze-out temperature = (1 + x~ A ) l / & Tw is larger than Tw, 
whereas the time scales as t' w = tw/(^ + x -4 ) 5 / 6 . In addition, 7/ is different from 7/ ~ 5 x 10~ 10 . 
However, since T/v depends on baryon density only logarithmically (see (|20|)), the temperature 
T' N remains essentially the same as Tn, while the time t' N scales as t' N = /{I + x^ 4 ) 1 / 2 . Thus, 
for the mirror 4 He mass fraction we obtain: 

y/~oy/^ \ 2exp[-t jV /T(l+x- 4 ) 1 / 2 ] 

r 4 _ 2x n (t N ) - i - exp[Am/TV(1 - x _ 4)1/6] • (22) 

We see that I4 is an increasing function of In particular, for x — > one has Yl — > 1. 



In reality, eq. (22) is not valid for small x, since in this case deuterium production through 
reaction n + p <-> d + 7 can become ineffective. By a simple calculation one can make sure that 
for x < 0.3 • (7/ x lO 10 ) -1 / 2 , the rate at which neutrons are captured to form the deuterium 
nuclei, T' N = n' p aN ~ v' n '~/ a N, where ~ T' 3 is the M-photons density and — 4.5 • 10 -20 
cm 3 s _1 is the thermal averaged cross section, becomes smaller than the Hubble rate H(T') 
for temperatures T' > T' N . In this case M- nucleosynthesis is inhibited, because the neutron 
capture processes becomes ineffective before deuterium abundance grows enough to initiate the 
synthesis of the heavier elements. Therefore, for any given 7/, first increases with increasing 
1/x, reaches a maximum and then starts to decrease. The true dependence of Yl on the x 
computed for different values of 7/ by standard BBN code (20], is presented in Fig. 0. The 
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Figure 3: The primordial mirror 4 He mass fraction as a function of x. The dashed curve represents the 
approximate result of cq. (pl|). The solid curves obtained via exact numerical calculation correspond, 
from bottom to top, to rf varying from 1CP 10 to 10 -6 . 



Hubble expansion rate of the mirror world was implemented, for each value of x, by taking an 
effective number of extra neutrinos as AN„ = 6.14 • x~ 4 . 

We have to remark, however, that in the most interesting situation when (5 = Qg/Clg = 
x 3 rj' ' jr\ > 1, the condition x < 0.3 • (7/ x lO 10 ) -1 / 2 is never fulfilled and the behaviour of Y£ is 
well described by the approximate formula (p2|). Hence, in this case I4 is always bigger than 
I4. In other words, if dark matter of the universe is represented by the baryons of the mirror 
sector, it should contain considerably bigger fraction of primordial 4 He than the ordinary world. 



5 Mirror baryons as dark matter 

We have shown that mirror baryons could provide a significant contribution to the energy density 
of the universe and thus they could constitute a relevant component of dark matter. Immedi- 
ate question arises: how the mirror baryon dark matter (MBDM) behaves and what are the 
differences from the more familiar dark matter candidates as the cold dark matter (CDM), 
the hot dark matter (HDM) etc. In this section we briefly address the possible observational 
consequences of such a cosmological scenario. 

In the most general context, the present energy density contains relativistic (radiation) 
component Q r , non-relativistic (matter) component O m and the vacuum energy density £l\ 
(cosmological term). According the inflationary paradigm the universe should be almost flat, 
CIq = J7 m + Q r + $7 a ~ 1 5 which well agrees with the recent results on the CMB anisotropy [21]. 



The Hubble parameter is known to be Hq = lOOh km s 1 Mpc 1 with h = 0.6 — 0.8, and for 
redshifts of the cosmological relevance, 1 + z = T/Tq S> 1, it becomes 

H(z) =H [a r (l + z) A + n m (l + z) 3 \ . (23) 

In the context of our model, the relativistic fraction is represented by the ordinary and mirror 
photons and neutrinos, ft r h 2 = 4.2 x 10~ 5 (1 + x 4 ), and contribution of the mirror species is 
negligible in view of the BBN constraint x < 0.64. As for the non-relativistic component, it 
contains the O-baryon fraction and the M-baryon fraction £l' B = /3$7_b, while the other types 
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of dark matter, e.g. the CDM, could also present and so Q m = £Ib + ^'b + ^cdm-Q At present 
it is not completely clear what is the matter fraction in the universe. Many observational data 
favour Q m ~ 0.3 while the rest of the energy density is due to the cosmological term, ~ 0.7, 
but also Q m ~ 1 cannot be excluded. 

The important moments for the structure formation are related to the matter-radiation 
equality (MRE) epoch which occurs at the redshift (we denote (f2 m /i 2 )o.2 = (r2 m /i 2 /0.20)): 



1 + = ^ « 2.4 • io 4 i^ = 4800 x (^ 2 )o.2 (24) 



and to the plasma recombination and matter-radiation decoupling (MRD) epochs. The latter 
takes place only after the most of electrons and protons recombine into neutral hydrogen and 
the free electron number density n e diminishes, so that the photon scattering rate T 7 = n e aT = 
X e rjn^aT drops below the Hubble expansion rate H(T), where ot = 6.65 • 10 -25 cm 2 is the 
Thomson cross section. In condition of chemical equilibrium, the fractional ionization X e = 
n e /nB is given by the Saha equation, which for X e <C 1 reads: 

X^(l-Y*) 1/2 ^(^X VA e- B/2T (25) 



rj 1 / 2 \m 



'<(:' 



where B = 13.6 eV is the hydrogen binding energy. Thus we obtain the familiar result that in our 
universe the MRD takes place in the matter domination period, at the temperature T dec ~ 0.26 
eV which corresponds to redshift 1 + Zdec = 

T dec /T ~ 1100. 

The MRD temperature in the M-sector T dec can be calculated following the same lines as 
in the ordinary one. Due to the fact that in either case the photon decoupling occurs when the 
exponential factor in eq. (|2^) becomes very small, we have T dec ~ Td ec , up to small logarithmic 
corrections related to rj , different from r], Y4. Hence 

1 + z' dcc ~ x- l {l + z dcc ) ~ 1.1 • lO 3 ^ 1 (26) 

so that the MRD in the M-sector occurs earlier than in the ordinary one. Moreover, for x 
less than x eq = 0.23(O m /i 2 )^2, the mirror photons would decouple yet during the radiation 
dominated period (see Fig. ||). 

Let us discuss now cosmological evolution of the MBDM. The relevant length scale for the 
gravitational instabilities is characterized by the mirror Jeans scale X'j ~ v^ir/Gp) 1 ^ 2 , where 
p(z) is the matter density at a given redshift z and v' s {z) is the sound speed in the M-plasma. 
The latter contains more baryons and less photons than the ordinary one, p' B = f3ps and 
p'^ = x 4 py. Let us consider for simplicity the case when dark matter of the universe is entirely 
due to M-baryons, Q m ~ £l' B . Then we have: 



c (1 , M 



-1/2 



v' s (z) ~ — 1 + 



B 



1 A P'*y I \/3 



3 A 1 \ 1 + z, 



1 + + ^ l + z 



cq 



-1/2 



(27) 



Hence, for the redshifts of the cosmological relevance, z ~ z e q, we have v' s ~ 2x 2 c/3 < c/y/3, 
quite in contrast with the ordinary world, where v s ~ c/ v3 practically till the photon decoupling, 
z = 1100. 

The M-baryon Jeans mass M'j = ^p m \'j reaches the maximal value at z = z' dcc ~ 1100/x, 
M'j{z' dec ) ~ 6-10 17 xx 6 [l+(a; e q/x)]- 3 / 2 (Q m /i 2 )^^ M . Notice, however, that Mj becomes smaller 

6 In the context of supersymmetry, the CDM component could exist in the form of the lightest supersymmetric 
particle (LSP). It is interesting to remark that the mass fractions of the ordinary and mirror LSP are related as 
f^Lgp — rrilLSP- In addition, a significant HDM component Sl„ could be due to neutrinos with order eV mass. 
The contribution of the mirror neutrinos scales as = x 3 Q v and thus it is irrelevant. 



11 



10 3 



10" 2 10" 1 1 

x 

Figure 4: The M-photon decoupling redshift 1 + z' dec as a function of x (solid). The long-dash line marks 
the ordinary decoupling redshift 1 + z<\ C c = 1100. We also show the matter-radiation equality redshift 
1 + z oq and the mirror Jeans-horizon mass equality redshift 1 + z' c , for the cases ft m h 2 = 0.2 (respectively 
lower dash and lower dott) and fl m h 2 = 0.6 (upper dash and dott). 

than the Hubble horizon mass Mh = ^pH~ s starting from a redshift z c = 750x~ 4 (f2 m /i 2 )o.2) 
which is about z eq for x = 0.64, but it sharply increases for smaller values of x (see Fig. 4). 
So, the density perturbation scales which enter horizon at z ~ z eq have mass larger than M'j 
and thus undergo uninterrupted linear growth immediately after t = t eq . The smaller scales for 
which M'j > Mh instead would first oscillate. Therefore, the large scale structure formation is 
not delayed even if the mirror MRD epoch did not occur yet, i.e. even if x > x cq . The density 
fluctuations start to grow in the M-matter and the visible baryons are involved later, after being 
recombined, when they rewritte the spectrum of already developed mirror structures. 

The main feature of the MBDM scenario is that the M-baryon density fluctuations should 
undergo the strong collisional damping around the time of M-recombination. The photon diffu- 
sion from the overdense to underdense regions induce a dragging of charged particles and wash 
out the perturbations at scales smaller than the mirror Silk scale X' s ~ 10 x /(x)(f2 m /i 2 ) ?/ 4 

Mpc, where f(x) = x 5//4 for x > x eq , and f(x) = (x/x eq ) 3//2 Xeq 4 for x < x eq . 

Thus, the density perturbation scales which can run the linear growth after the MRE epoch 
are limited by the length A^. This could help in avoiding the excess of small scales (of few Mpc) 
in the large scale power spectrum without tilting the spectral index.]] The smallest perturbations 
that survive the Silk damping will have the mass M' s ~ x 15 ' 4 (n m h 2 ) 2 /4 10 13 M , which should 
be less than 2 x 10 12 M Q in view of the BBN bound x < 0.64. Interestingly, for x ~ x eq we have 
M' s ~ 4 x 10 10 (ri m /i 2 )Q2 M , a typical galaxy mass. 

To some extend, the cutoff effect is analogous to the free streaming damping in the case of 
warm dark matter (WDM), but there are important differences. The point is that alike usual 
baryons, the MBDM should show acoustic oscillations whith an impact on the large scale power 
spectrum. In particular, it is tempting to imagine that the M-baryon oscillation effects are 
related to the anomalous features observed at 100/i _1 Mpc clustering p^ j. 

In addition, the MBDM oscillations transmitted via gravity to the ordinary baryons, could 
cause observable anomalies in the CMB angular power spectrum for Vs larger than 200. This 
effect can be observed only if the M-baryon Jeans scale Aj is larger than the Silk scale of ordinary 
baryons, Ag ~ 3(3 1 ^ 2 (Q m h 2 ) ^ 4 Mpc, which sets a principal cutoff for CMB oscillations around 
I ~ 1200. As we have seen above, this would require enough large values of x, near the upper 

7 One should keep in mind also delaying of the MRE moment, e.g. by taking Q. m h 2 ~ 0.2. This would 
correspond to the case when M-baryons provide only a fraction of the present energy density, say Q' B ~ 0.3, while 
the rest is due to the cosmological term, S7a ~ 0.7. 
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bound set by the BBN constraints: x ~ 0.6 or so. The detailed analysis of this effect will 
be given elsewhere. In our opinion, together with the possible effects on the large scale power 
spectrum, it can provide a most direct test for the MBDM and can be verified by the next CMB 
experiments with higher sensitivity. 

Clearly, for small x the M-matter recombines before the MRE moment, and thus it should 
rather manifest as the CDM as far as the large scale structure is concerned. However, there 
still can be crucial difference at smaller scales which already went non-linear, like galaxies. 
In our scenario, dark matter in galaxies and clusters can contain the mixed CDM and MBDM 
components, and can be even constituted entirely by the mirror baryons. Then one can question 
whether the MBDM distribution in halos can be different from that of the CDM? Namely, 
simulations show that the CDM forms triaxial halos with a density profile too clumped towards 
the center, and overproduce the small substructures within the halo. As for the MBDM, it 
constitutes a sort of collisional dark matter and thus potentially could avoide these problems, 
at least clearly the one related with the excess of small substructures. 

The halo distribution in galaxies depends on the mass M and on the self-scattering cross 
section a of dark matter. In our case it mainly consists of the mirror hydrogen atoms, and so 
Mh — 1 GeV and an ~ 10 -16 cm 2 . At the first glance, this is in strong discrepancy with the 
range a/M ~ 10~ 23 — 10~ 24 cm 2 /GeV preferred by the analysis of ref. 

However, one has to take into account the possibility that during the galaxy evolution the 
bulk of the M-baryons could fastly fragment into the stars. A difficult question to address 
here is related to the star formation in the M-sector, also taking into account that its tempera- 
ture/density conditions and chemical contents are much different from the ordinary ones. In any 
case, the fast star formation would extint the mirror gas and thus could avoide the M-baryons to 
form disk galaxies as ordinary baryons do. The M-protogalaxy, which at certain moment before 
disk formation essentially becomes the collisionless system of the mirror stars, could maintain 
a typical elliptical structure.^ Certainly, in this consideration also the galaxy merging process 
should be taken into account. As for the O-matter, within the dark M-matter halo it should 
typically show up as an observable elliptic or spiral galaxy, but some anomalous cases can be 
also possible, like certain types of irregular galaxies or even dark galaxies dominantly made out 
of M-baryons. 

The central part of halo can nevertheless contain a large amount of the ionized mirror gas 
and it is not excluded that it can have a quasi-spherical form. Even if the stellar formation is 
very efficient, the massive mirror stars in the dense central region fastly evolve^ and explode as 
supernovae, leaving behind the compact objects like neutron stars or black holes,|^] and repro- 
ducing the mirror gas and dust. Although the cross section ajj is large, it does not necessarily 
implies that the galaxy core will collapse within a dynamical time, since the inner halo should 
be opaque for M-particles. They undergo many scatterings and escape from the system via 
diffusion, so the energy drain can be small enough and the instability time can substantially 
exceed the age of the universe |2jJ. 

In the galactic halo (provided that it is the elliptical mirror galaxy) the mirror stars should 

8 In the case of asymmetric M-world with ( = i//d > 1, the mirror electron mass should scale as m' e ~ £ro e 
while the nucleon mass remains ~ 1 GeV. Therefore, the elastic scattering cross-section of the M-hydrogen atoms 
scales as a' H ~ C~ 2cr B and so for £ 3> 1 one could obtain a reasonably small cross section. 

9 In other words, we speculate on the possibility that the M-baryons form mainly the elliptical galaxies. For 
a comparison, in the ordinary world the observed bright galaxies are mainly spiral while the elliptical galaxies 
account about 20 % of them. Remarkably, the latter contain old stars, very little dust and show no sign of active 
star formation. 

10 Since mirror matter contains more helium, the mirror stars should evolve much faster than the ordinary ones. 

11 Another tempting issue is whether the M-matter itself could help in producing big central black holes, with 
masses ~ 10 7 M Q , which are thought to be main engines of the active galactic nuclei. 
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be observed as Machos in gravitational microlensing ||, ||]. Leaving aside the difficult question 
of the initial stellar mass function, one can remark that once the mirror stars could be very old 
and evolve faster than the ordinary ones, it is suggestive to think that most of massive ones, 
with mass above the Chandrasekhar limit Mq^ ~ 1.5 M & have already ended up as supernovae, 
so that only the lighter ones remain as the microlensing objects.^ The recent data indicate the 
average mass of Machos around M ~ 0.5 M & , which is difficult to explain in terms of the brown 
dwarves with masses below the hydrogen ignition limit M < 0.1M Q or other baryonic objects 
]. Perhaps, this is the first observational evidence of the mirror matter? 
It is also plausible that in the galactic halo some fraction of mirror stars exists in the form 
of compact substructures like globular or open clusters. In this case, for a significant statistics, 
one could observe interesting time and angular correlations between the microlensing events. 

6 Conclusions 

We have discussed cosmological implications of the parallel mirror world with the same micro- 
physics as the ordinary one, but having smaller temperature, T 1 < T, with the limit T'/T < 0.64 
set by the BBN constraints. Therefore, the M-sector contains less relativistic matter (photons 
and neutrinos) than the O-sector, £l' r <C Sl r - On the other hand, in the context of the GUT or 
electroweak baryogenesis scenarios, the condition T < T yields that the mirror sector should 
produce a larger baryon asymmetry than the observable one, rj' B > ?]bJ3 So, in the relativis- 
tic expansion epoch the cosmological energy density is dominated by the ordinary component, 
while the mirror one gives a negligible contribution. However, for the non-relativistic epoch the 
complementary situation can occur when the mirror baryon matter density is bigger than the 
ordinary one, £l' B > Q,b- Hence, the MBDM can contribute the dark matter of the universe 
along with the CDM or even constitute a dominant dark matter component. 

We have also shown that the BBN epoch in the mirror world proceeds differently from 
the ordinary one, and it predicts the mirror helium abundance in the range = 0.4 — 0.8, 
considerably larger than the observable Y4 ~ 0.24. 

Unfortunately, we cannot exchange the information with the mirror physicists and combine 
our observations. (Afterall, since two worlds have the same microphysics, the life should be 
possible also in the mirror sector.) However, there can be many possibilities to disentangle the 
cosmological scenario of two parallel worlds with the the future high precision data concerning the 
large scale structure, CMB anisotropy, structure of the galaxy halos, gravitational microlensing, 
oscillation of the neutrinos or other neutral particles into their mirror partners, etc. 
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